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Certain spacelike mesonic correlation lengths serve as interesting theoretical probes 
for the reliability of perturbation theory in high-temperature QCD, are directly 
sensitive to chiral symmetry restoration and to the axial anomaly, and might also 
have indirect phenomenological signatures in idealised heavy ion collision experi- 
ments. I review here the weak-coupling predictions for some of these correlation 
lengths, to be compared with results from lattice Monte Carlo simulations. 
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1. Introduction 

<D ! 

e-j , At a finite temperature T, Lorentz symmetry is broken by the heat bath, so 

that temporal and spatial directions are, in general, unrelated. Correlation 
functions in the (Minkowski) time direction can be used to define, after a 
Fourier transform, spectral functions for various operators, and the spec- 
■ tral functions in turn determine fundamental "real-time" properties of the 

plasma, such as particle production rates. 1 These are then directly measur- 
able in heavy ion collision experiments (modulo problems with background 
contamination, etc). 

Correlation functions in spatial directions, in contrast, address ques- 
tions such as: At which length scales are thermal fluctuations correlated? 
At which length scales are external charges screened? In principle these 
"static" observables are also physical and may lead to detectable signals, 
however in practice the relations are indirect and therefore weaker than 
for real-time observables. On the other hand, static quantities are emi- 
nently suited to measurements in lattice experiments. Lattice simulations 
will hence be the most immediate "phenomenological" reference point in 
the following. 

There is naturally a vast variety of different operators that can be cor- 
related. The operators can be classified according to their discrete and 
continuous global symmetry properties, leading to many independent cor- 
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relation lengths. The analytic structures of various purely gluonic correla- 
tors were discussed in Ref. 2, and the corresponding lattice measurements 
are rather precise by now. 3-8 The purpose of this talk is to analyse the 
structures of mesonic correlators, for which the status has been somewhat 
less advanced. This talk is based on the original study in Ref. 9. 

2. Detailed setup 

As just mentioned, we will focus on correlation lengths £ of mesonic ob- 
servables, that is operators of the type O = $ TF a ^, where 

r = {1,75,7^.7^75} , (1) 
pa = { F ^F n }, F s = I iVf x Nf , F n = traceless N f x N f matrix , (2) 

iVf is the number of quark flavours, and all colour, spinor and flavour in- 
dices have been suppressed. We may recall that at least some among these 
operators have direct physical significance, for instance 

V'75-F lS ' ! /' oc ry'-meson , (3) 

■ipy 5 F n ip oc pion , (4) 

■tpjoF s ip oc baryon number density , (5) 

vp7oF n ip oc electric charge density (for Nf = 3) . (6) 

The lattice determination of spatial correlation functions for operators of 
this type was pioneered by DeTar and Kogut 10 and others 11 a long time 
ago, but quantitatively significant results relevant for the physical infinite 
volume chiral continuum limit, are only being produced presently. 12-15 In 
order to streamline the discussion, we will restrict to one of the operators 
in the following, namely the interpolating field for the pion: ir ~ i'ijjj^F n tp. 
For ^7o-F s V', see Sec. 4, and for the other operators, Ref. 9. 

The spatial correlation function of the pion field has at large distances 
(|x| — > oo) the structure 

C x = £ /T dT(n(T,xM0,O))~^exp(--^ , (7) 

where r is the Euclidean time coordinate of the imaginary time formulation. 
In Eq. (7), n is some unspecified power, while £(T) defines the correlation 
length we are interested in. We will denote the inverse of the correlation 
length by m = £ _1 , and call it a "screening mass". 

The question now is, how does m(T) behave as the temperature is in- 
creased from below to above the critical temperature T c of the chiral phase 
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transition? At small temperatures, chiral symmetry is broken, and the pion 
is massless in the chiral limit, resulting in an infinite correlation length, or 
to = 0. At high temperatures, on the other hand, the screening mass ap- 
proaches 2nT, as we will recall presently. Thus the pion screening mass 
provides for a finite and gauge-invariant "order parameter" for chiral sym- 
metry restoration. 

The reason that the screening mass equals 2nT at high temperatures, 
is easily understood. For T ^> T c , asymptotic freedom sets in, and the 
correlator can be determined in perturbation theory. A computation of the 
leading order graph, 




and a subsequent Fourier transform [C p = J exp(— iu> n T— ip-x)C x ], yield 

n— — oo ' 1 

Now, the singularity closest to the origin in momentum space (determining 
the behaviour at the largest distances in configuration space) is seen to be 
a two-particle threshold at |p| = ±2iujQ. Referring for a moment to the 
spatial directions in a (2+l)-dimensional language, with X3 as the "time" 
coordinate, both quarks are on-shell at this point and have a minimal "en- 
ergy" ip3 = ±wq [the minimal energy appears for a momentum = 0, 
where we denoted p = (p±,p3)]. Thus, taking an inverse Fourier transform 
back to configuration space, we obtain the exponential decay advertised 
above, C x ~ exp(— 2wo|x|). 

3. Next-to-leading order for flavour non-singlet correlators 

The question we would like to address next is, what is the first weak- 
coupling correction to to = 2nT? This computation has been described in 
Ref. 9, and we only reiterate the main steps here. 

For flavour non-singlet operators, the relevant graphs might appear to be 
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but in fact there are infinitely many higher order graphs that need to be 
taken into account, symbolically of the types 




These are not suppressed, because near the two-particle threshold, the 
quarks are almost on-shell, p 2 = u> 2 + + p 2 « 0, and consequently 
the dimensionless expansion parameter associated with adding a further 
line can be seen to be 0(g 2 T/\]& |) ~ 0(1). In other words, the situation is 
analogous to the computation of the energy of a two-quark bound state at 
zero temperature, where again infinitely many graphs contribute. 

The analogy with the bound-state system suggests immediately a tool 
for organising the computation: As in the study of quarkonia at T = 0, 
energies of bound states can be addressed with an effective theory called 
"Non-Relativistic QCD" (NRQCD). 16 NRQCD techniques were first em- 
ployed in the present context by Huang and Lissia. 17 

In order for an effective description to apply, the system must possess a 
scale hierarchy. Let us recall why one exists here, at high enough tempera- 
tures where the QCD coupling constant g is small. The basic point is that 
because quarks of a definite Matsubara frequency w n interact with bosonic 
Matsubara zero-modes only, we expect that their off-shellness is related to 
the momentum scales of the latter: \ip% ± uj n \ < gT. More precisely, we may 
recall that the static potential has in (2+1) dimensions the structure 

F(x ± )~ 5 2 Tln|xJ . (9) 

Thus the typical (transverse) momentum p± of the bound-state con- 
stituents satisfies 

p 2 Ju n ~ V(x x ) => pi ~ (9T) 2 , (10) 
and the binding energy is (dropping again possible logarithms) 

\i P3 ±U n \~p 2 ± /u n ~V(x ± )~g 2 T. (11) 

We thus indeed find a scale hierarchy 

\ip 3 ±Lu n \ < |p x | <w„ ■ (12) 

The effective action for this kinematic range can most easily be con- 
structed by choosing a convenient basis for Dirac matrices (making 7073 
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diagonal), and rewriting then the Dirac spinor as 

d3) 

where \, <j> are two-component spinors. In this basis the pion operator 
becomes ir ~ i(x — 0'(73x)j where 03 is a Pauli matrix. Restricting to 
the Matsubara mode luq and expanding the QCD action to order 1/luq, the 
(on-shell) effective Lagrangian reads 

C = i*t (uo - gA + D 3 - x + itf (uo ~ gA -D 3 - , (14) 
where D3 is a covariant derivative. We may note that: 

• Both Aq and A3 play an important dynamical role. 

• The transverse gauge fields A\ , A2 , on the other hand, can be ig- 
nored, as long as we are interested in an energy shift of 0(g 2 T): 
they are of higher order than Vj_ ~ gT. (They would be of order 
unity with respect to the magnetic scales Vj_ ~ g 2 T, but that gives 
an energy correction of order 0(g A T) only.) 

• To be consistent at 0(g 2 T), we should replace ojo of the tree-level 
effective Lagrangian by a matching coefficient M = ojq + 0(g 2 T) 
that needs to be determined. 

While the final value of the matching coefficient M is unambiguous, 
there are many ways to determine it. In order to avoid computing wave 
function normalisation factors, we determine M by matching on-shcll quark 
self-energies in the original QCD as well as in the effective theory of Eq. (14). 
On the side of QCD, a simple 1-loop graph, 



p 2 = 

produces near the pole the inverse propagator 

7m -i)nn 



E(p) = # - ig 2 C F 



(p - q) 2 (q 2 + X 2 ) b 



(15) 



where Cp = (N 2 — 1)/2N C , N c is the number of colours, S is the standard 
imaginary-time integration measure, (...)/ and (...)h denote fermionic and 
bosonic Matsubara four-momenta, respectively, and we have introduced a 
gluon mass A as an intermediate infrared (IR) regulator. It can be seen 



February 2, 2008 0:49 Proceedings Trim Size: 9in x 6in 



sewm04 



6 



that the integral in Eq. (15) remains finite for A — > 0, and then, up to order 
0(g 2 ) 1 solving for the zero of Eq. (15) is equivalent to solving 



0, 



(16) 



where m 2 s = j rCp/l is nothing but the "hard effective mass" of quarks 
introduced recently in Ref. 18. 

On the side of NRQCD, the computation of the quark self-energy is to 
be carried out order by order in l/w . 19 It is easy to see that, in fact, 1-loop 
corrections then vanish, a such that the on-shell point is determined directly 
by the mass scale M appearing in the tree- level propagator, ip 3 = ±M. 

Consequently, solving for the zeros in ip 3 of Eq. (16) to order 0(g 2 T), 
we see that M needs to match oj + m c 2 ff /2w > leading to 



M = lu + g'T^- + 0(g 4 T) . 



(17) 



This value replaces uj in the first terms inside the parentheses in Eq. (14). 

The parameters of the effective theory having been determined, it re- 
mains to solve its dynamics to order g 2 T. We define the correlation function 

C(r, x 3 ) = J (f (R + |, z 3 ) <?3X (R - \ x 3 ) X f (0)^(0)) , (18) 

and integrate out A ,A 3 . To order 0(g 2 ), C(r,x 3 ) satisfies the partial 
differential equation 

1 



where 



V(r) 



+ 2M V; 

LU 



+ V(r) C(r,x 3 ) ot6(x 3 )5(r) , 



fTC F 



f d 2 - 2e g fl 
mar 



= g 2 T 



Cf_ 
2ir 



In 



- e l * r 1 + e iC i 
q 2 + X 2 q 2 + m 



D ) 



+ -f E - K a (m D r) 



(19) 



(20) 
(21) 



On the last line we took the continuum limit (e — > 0) and removed the IR 
regulator from the propagator of A 3 (A — ► 0). Moreover, Kq is a modified 
Bcssel function, and is the Debye mass appearing in the propagator of 
A , m 2 D = g 2 T 2 (N c /2> + Nf/6). Solving for the lowest energy level of the 
Schrodinger equation following from Eq. (19) at x 3 ^ 0, we finally obtain 



2nT 



g 2 T 



1 



cw 

2?r V2 



■•) 



^0 



0.164 



(22) 



a It is important to note that at this point we are using un-resummed gluon propagators 
both on the QCD and on the effective theory side. 
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Figure 1. The pion screening mass as a function of the temperature. The dashed line is 
a guess for how the low-T and high-T limits might be connected, and in the continuous 
form drawn assumes a second order transition (i.e. Nf = 2). Note, however, that no sign 
of the indicated overshooting of 27rT has been observed in lattice simulations. 



where Eq depends mildly on Nf (cf. Ref. 9), and we have shown explicitly 
the value relevant for the quenched theory. b 

Note that in Eq. (22), the factor 1/2 inside the parentheses comes from 
the "constituent mass" correction in Eq. (17), while E comes from the 
solution of the Schrodingcr equation. If the solution corresponds to a bound 
state, why then is E positive? The answer lies simply in the form of the 
potential in Eq. (21): it is a logarithmically "confining" potential, with 
V(r — > 0) — > — oo, V(r — > oo) — » +oo, whose overall normalisation just 
happens to be such that the lowest energy eigenvalue is positive. 

The result of Eq. (22) is plotted, with a certain choice of renormalisation 
scale for g 2 (cf. Ref. 3), in Fig. 1. We note that the correction appears small 
even at realistic T, but positive in sign. The change of m across T c is thus 
very rapid — it is not smoothed by the 0(g 2 ) term like for the pressure, 20 
but rather made more pronounced. Close enough to T c the mass of course 
has to decrease, either continuously (the case for Nf — 2, indicated with a 
dashed line in Fig. 1), or with a small discontinuity (for Nf — 0, 3). 

For the pressure, the next-to-next-to- leading order correction, 0(g 3 ), 
comes with an opposite sign to 0(g 2 ), and appears very large. 21 (At the 
same time the ultimate result, including higher orders still, 22 resembles 



b To be precise, the quenched theory value applies in the "perturbative" broken Z(N C ) 
vacuum, i.e. the one where the phase of the Polyakov loop is trivial. In the unquenched 
theory this vacuum constitutes in any case the global ground state of the system. 



February 2, 2008 0:49 Proceedings Trim Size: 9in x 6in 



scwm04 



8 



more the 0(g 2 ) approximation.) It would be interesting to find out whether 
there is an 0(g 3 ) correction in the pion correlation length as well. 

4. Flavour-singlet correlators 

For flavour singlets the pattern is very different from the flavour non- 
singlets, as the singlet operators couple to purely gluonic "glueball" states. 
For baryon density, for instance, one finds that 8,23 



or graphically, that the correlation function is dominated by graphs like 



The corresponding screening mass has been determined some time ago, 
and evaluates to m(Tr[Ag]) w 5T at T « 2T C . The inverse, £ « 1/5T, is 
then the distance scale at which fluctuations in baryon number density are 
correlated. The glueball operators mediating the correlations of a number 
of other mesonic flavour-singlets have been worked out in Ref. 9. 

5. On Ua(1) symmetry at high temperatures 

The U,a(1) axial symmetry is broken explicitly by the anomaly. It has 
been suggested, 24 however, that it might get "effectively restored" some- 
what above the critical temperature, say at T>2T C . What is meant by 
this is that the topological susceptibility, x = { l/2 )/V, where v is the topo- 
logical charge of the gauge configuration, or the net number of zero-modes 
of the Dirac operator, and V is the volume, decreases rapidly at tempera- 
tures above T c . This becomes quite obvious particularly if the number of 
colours is increased above N c = 3; see, e.g., Ref. 25 and references therein. 
Consequently, it could be expected that anomalous effects cease to operate. 

It may be appropriate to remark at this point, though, that it is possi- 
ble to find observables for which the Ua(1) symmetry appears not to get 
restored. Consider, for instance, the screening masses for the operators 
V$ = V'73'0j ^3 — V'7375'0! measured from correlation functions in the 
X3-direction. Now, Ux(l) symmetry restoration would imply that the cor- 
relation functions for V£, A$ are identical and, therefore, that the screening 
masses are identical. Yet, a perturbative analysis can be used to indicate 9 




(23) 
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that A§ couples to the glueball-like state g 2 tij3 Tr[A Fy] of the dimen- 
sionally reduced theory, 26 and has thus a non-vanishing screening mass 
m = mu + 0(g 2 T), 2 ' 3 where the next-to-leading order correction is also 
known numerically. 5 At the same time, the correlator related to V 3 S is ex- 
actly conserved, and has no non-trivial screening mass associated with it. 
Therefore, the Ua(1) symmetry appears broken even at high temperatures. 

6. Conclusions 

Various mesonic correlation lengths are well-defined gauge-invariant phys- 
ical observables that seem to be computable, for flavour non-singlets, at 
least up to next-to-leading order in the gauge coupling g 2 , as outlined in 
this talk. Higher orders could also be reached in principle. For flavour 
singlets, on the other hand, the correlation lengths can be related to those 
of purely gluonic states, which have been determined previously from the 
dimensionally reduced theory with good numerical precision. 4 ' 5 ' 8 

The comparison of these analytic predictions with results from four- 
dimensional lattice simulations poses an interesting puzzle. Indeed, lattice 
results so far are consistent with the leading order value m = 2nT, as soon 
as the temperature is a bit above the critical one, 12-15 while the next-to- 
lcading order correction advertised above is positive, requiring an "over- 
shooting" of m = 2nT. Lattice results involve systematic errors, though, 
so that the discrepancy cannot be considered too serious at present. 

On the side of analytic efforts, one could envisage a number of compu- 
tations which might facilitate the comparison with ever more precise future 
simulations. For instance, the leading discretization effects could be de- 
termined analytically, even though this of course is specific to the fermion 
action used. (Different groups have rather different preferences here.) An- 
other point of possible relevance is that some of the correlation functions 
involve an ^-dependent prefactor in front of the exponential, and it might 
be useful to take it into account in the fitting procedure. 

To conclude, let us recall that the ultimate theoretical goal of the ex- 
ercise discussed here is to estimate the reliability of perturbation theory 
at finite T, in order to then use perturbation theory with more confidence 
for other observables for which lattice is not well suited. Given that good 
QCD predictions can eventually be obtained for static correlation lengths 
and that there is a remarkably stark difference between the two phases, 
one is however also lead to wonder, once again, whether some interesting 
phenomenological signatures might be found for these quantities. 
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